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Abstract 

We establish the relation between full traces of the Green functions for 
some initial and the Darboux transformed one-dimensional two component 
Dirac problems with the most general form of potential. The result is used 
to check the completeness of set of wave functions obtained by the Darboux 
transformation of the eigenfunctions set for the initial Dirac problem with 
some typical boundary conditions. 
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1 Introduction 

Supersymmetric quantum mechanics (SQM) [HIS] provides an interesting framework 
within which to analyse quantum problems. In particular, it allows one to investi- 
gate the spectral properties of a wide class of quantum models and to generate new 
systems with given spectra. SQM gives new insight into the problem of spectral 
equivalence of Hamiltonians, which, historically, has been constructed as a factor- 
ization method in quantum mechanics [3] and as Darboux-Crum transformations in 
mathematical physics f^. 

The Darboux transformation method for the one- dimensional stationary Dirac 
equation is equivalent to the underlying quadratic supersymmetry and the factor- 
ization properties of the Dirac equation [5]. Application of this method to the Dirac 
equation is studied in the papers P, [7] . 

In the previous paper [7j we have established the connection between Green 
functions of initial and Darboux transformed one-dimensional two component Dirac 
equations for the case of especial matrix structure of the potential. Here we consider 
the same problem for the arbitrary matrix structure of the interaction Hamiltonian. 
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2 Green function of initial problem 

We consider the Dirac equation 

ho{x)tlj{x) = E^{X), V = {^l:^2f (1) 

with the Hamiltonian of the form 

ho = ta2d, + V{x), (2) 

Vo{x) = uj{x)I + (m + S{x))ai + q{x)ai, (3) 

where u;(x), S{x) and g(a;) are real functions of x\ m is the mass of a particle; di, (T2, 
(Ts are usual Pauli matrices. This equation have two linearly independent solutions. 

Denote them by ip and ^p. Introduce the Wronskian of this solutions with the 
help of equation 

W{i){x),ip{x)} = ^i{x)lp2{x) - ^2{x)ipi{x). (4) 
Let us prove that it doesn't depend on x. For this aim represent W in the form: 



W ^ ij 'yep, 7 = i(T2 = _^ Q J , 7 = -7- (5) 
Then 

w = i^^y-f^ + ^^j^' = tp^j^' -i-fij'ff- (6) 

Taking into account the Dirac equation 

7(^' = Eip-V^, (7) 
7^' = Eip-Vip, (8) 

we have 

W - ij'^{Eip-Vip)-{Eilj'^ -{Vijf)(p = ij'^{E-V-E + V^)ip. (9) 
Since = V,W' = 0. 

Now we are in position to construct Green function of problem under consider- 
ation. It is the solution of inhomogeneous equation 

{H-E)G{x,y) = 6{x-y). (10) 

It is easy to check that 

^ V{x)'4}'^ {y)e{x -y)+ i^{x)ip'^{y)e{y - x) ^^^^ 

Besides this for G{x, y) the spectral representation 

G(x,,) = E„'^^M») (12, 

is valid. Here y?„(x) is complete orthonormal system of eigenfunctions of the Hq with 
eigenvalues En. 
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3 Darboux transformation of the Green function 

Let us to construct the 2x2 matrix u, consists of two two- component eigenfunctions 
V^m, fn2 of the Dirac Hamiltonian ho corresponding eigenvahies En^: 

u = (v9n^,v?„J. (13) 

Thus, the function -u is a solution of the Dirac equation 

Hqu = uA, A = diag{Enj^, En^). (14) 

The operator 

L = u^ = ^ (15) 

ax ax 

allows us to generate the solutions ipn of the transformed Dirac equation 

hiipn = Enipn, (16) 

from solutions of the initial Dirac equation. The Hamiltonian and corresponding 
solutions of the transformed Dirac equation are the following 

(^n = Lipn (17) 

hi = 7^ + ^i, Vi = Vo + [-i,u^u-\ (18) 

where (pn are solutions of the initial equation. 

It is easy to check that Lipm = 0, Lipn2 = 0. The spectrum of hi doesn't contain 
eigenvalue En^, En2- 

The question arises: are the eigenfunctions ipn {n ^ 1,2) form the full system. 
In order to answer to this question it is necessary to construct Green function Gi 

{hi-E)Gi = 6{x-y) (19) 

and to calculate the expression 

f 11 

A = tr nGi{x,x)-Go{x,x))dx + - - + - -. (20) 

If A = 0, the set of if) form the complete system. 

Early expression, similar to (120]) . was considered for Schrodinger equation in [S]. 
Evidently that 

Q^l^^ ^ ^{x)ip^{y)<S>{x - y) +^ij{x)^'^{y)Q{y - x) ^^i) 

W 

W = i'^{x)-f<^{x). (22) 
It can be proved (see Appendix A) that W = {E — EnJ{E — En2)W. 
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Consider 

tr{tp{x)^^{x)} = i'^{x)^{x) = {L^{x)fLip{x). (23) 
It is evidently that 

^ = Lil) = il)\x) - u^u~^i){x) = (24) 

and 

^^(x) = [(^-v)r«^, 

^^(x)<^(a;) = {u'^i;f)'u^(p = {i;'^{u-^fu''^Y-{i;^{u-^f{u^^)']. (25) 

Since u is a real matrix (this follows from reality of all coefficients of the Dirac 
equation rewritten in component form) = . The last form in the equation can 
be represent in the form 

V'^L+c^, L+ = -d^ - {u^u-^)+. (26) 

Taking into consideration that (p = Lip, factorization relation [5] 

L^L = {ho — Eni){hQ — -EnJ, 



and 

we have 
Thus, 



ho(f = Elf, 

^""(p = ,lj^(p + {E-Er,.,){E-E^,)^'^ip. (27) 



t.G.ix^x)-'^^^ i^^'fy + iE-E^,)iE-E^,){Vv>) 



W {E - E^,){E - E^,)W 

= '^£^ + tr{G{x,x)}. (28) 
W 

Then the trace of difference of the Green functions is as follows 

tr{Gi{x,x)-Go{x,x)\ = ^^J:' . (29) 

W 



Evidently that fl29l) can also be presented in the form 



Thus, 



tr[Gi{x,x)-Go{x,x)] = . (30) 

W 



tr{ip^<p-ilj^<p} = C, G' = 0. (31) 
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Let us prove that 

C={2E-En,-En,)W. (32) 
Taking into account the Dirac equations for if)^ and u 

^^P' = {E-V), -lip' = {E-V)ip, ^u^u'' = ^l-V, (33) 
where VL = uAu~^ it is easy to see that 

ip'^ip-^'^ip = ip'^K^f - E)-f'^ E)] 

= 2E^^7(^-?/'^(n^7 + 7fi)(^. (34) 

fi^7 + 7^ = -ftrQ 
trQ = En, + En, 
_ ^ ^2E - E„, - E„.,)W 

is proved. 

Taking into account this relation one can obtain 
tr / [Gi{x,x) — Go{x,x)]dx 

J a 



Since 
and 

the relation 



1 

+ 



E — En, E — En2 



+ ij'{bMb)-^'{aMa) (35) 



or 



r'' 11 

A = tr [Gi{x,x)-Goix,x)]dx + - + - 

Ja ^ni - ^ ^n2 - ^ 

= ^^(6)^(6) -^^(a)^(a). (36) 

If right side of this relation equal to zero the wave functions form complete system 
of eigenfunctions of hi. In opposite case this system is incomplete. 



4 Examples 

In this Section we consider some examples of the application of fl36l) to the Dirac 

problem with homogeneous boundary conditions on the finite interval x G [a,b]. 
More explicitly, we look for the solutions of the Dirac equation 

hip = Eip, (37) 

that satisfy the boundary conditions of type 

'?/'i(a) cosa + ^/'2(a) sina = 0, (38) 

^i(6)cos/5 + (^2(&)sin/3 = 0. (39) 
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We start with the consideration of the case of free initial Hamiltonian 

ho = ia2-^+mas (40) 

and consider four types of boundary conditions: (i) a = /3 = 0; (ii) a = P = 7i/2; 
(iii) a = 0, /3 = 7r/2; (iv) a = n/2, P^O. 

Without the loss of generahty we can put a — 0, b — —1. 

Case (i). The eigenfunction of Hq consists of two branches: (a) positive eigen- 
values and (b) negative eigenvalues. Positive eigenvalues have a form: 



E„ = ^/kl + m'^, kn^nn, n = l,2, ... (41) 
The corresponding eigenfunctions are 

(p\' = sin(A;„j;), (p^^ ^r— ■ (42) 

Negative eigenvalues are as follows: 

E = -^/kl + m'^, kn^nn, 0,1,2,... (43) 

Corresponding eigenfunctions are 

T(n) . /, \ , (n) kn COS [knX) , . 

= sm{knx), -02 , n = l,2,... (44) 

£/ + r/i 

(45) 

= 0, V^f = c, 

where C is an arbitrary constant. 

The independent solution of the equation 

ho^ = EjP, E^{E^,E^) (46) 

are 

V'l = shiikx), ^^^_ kco<kx) ^ k^VE^- m?, (47) 

E + m 

kcosikx — x) 

(fi = sin(te-fc), (^2 = ^T"^ ■ (48) 

H/ -\- Tfl 

Then the full trace of the Green function is 

/E 
Go(x,x)dx = — cot(A;)-^. (49) 
k k^ 

Taking into account the following relation (see Appendix B) 

I ^ 2k 

n=l 
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we get 

m \ — - 1 

■com-- = 2i^E^ 

n=l 
oo ^ 

^E'^-El ' E + m 

n=l 
oo / 

E 



E m 1 E — m 

k ^ ' kP' ^ k'^ — n^vr^ k'^ 

n=l 

1 

n=l 

1 



n=l 



E- En E- En E + m 



that is just the spectral representation of the tr Gq{x, x)dx. 

Now we consider the Darboux transformation. The solutions of transformed 
equation are following 

^ = -0' - UxU~'^iIj, (51) 

where u is the transformation matrix. For the construction of u we choose the pair 
of the eigenfunctions t/''-^'' and ip^^^ or pair ip^'^^ and ^p^^\ 
In the first case the transformation matrix reads 

_ / sin(A;ix) \ _^ _ / Aj^ cot /cio; 0\ , . 

U — \^ (J fci cos(fcia:) j , UxU ~ _ g J ■ I'^'^J 

In the second case it reads 

sin(A;ix) ^ „ „ -i _ cot kix 

1 cos(A;i 
Ei+m 

This leads to the following expressions for components of functions 

ip{x) = Lip, and <f{x) = Lip : 



^ = \ C -h^p^ )' =[{E,-m) )■ ^^2) 



xpi = k cos (kx) — IT cot{'n'x)sm{kx), (54) 

1P2 = {E-Ei)sm{kx) {or ip2 = {E - Ei))sm{kx), (55) 

cpi = k cos{kx — k) — TT cot{7cx) sin{kx — k), (56) 

(^2 = (^ - ^1) sin(A;x - A;) (or ^2 = (^ - ^1)) sin(A;x - A;). (57) 

It is interesting to note that at the left side {x = 0) of the interval both com- 
ponents of ip (V'l and 1P2) are zero, that is i^i^a) = 0. Similarly, <^i(l) = '^2(1) = 
or {p{b) = 0. Thus, the right side of eq. (IHUI) is zero and consequently system of 
function L(l)^^\ L0*^") is complete. 

Case (ii). The eigenspectrum again consists of two branches: 
(a) En = ^/kn + ra^i kn = mr, n = 0,1, ... 
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(b) En = -^/kl + w?, /cn = nvr, n = 0, 1,... 
Solutions '?/', if reads: 

Aj cosf Aj^I 

'4'2{x) = smikx), i'lix) = —- , k = V-E^ — rn^ (58) 

E — m 

LP2 = sm{kx-k), (fi{x) = ^^^^^^ (59) 

E — m 

Then the corresponding full trace of the Green function may be represent in the 
form: 



-. (60) 

m 



The transformation matrix u is 



/ 1 kiCos(kix)/(Ei — m) \ ( \ kiCOs(kix)/(Ei — m) \ 

M = „ ■ \ ] or u = { ■ \ (61) 

Y U sm^fcix) J y sm[kix) J 

and components of solutions of the transformed Dirac equation are: 

■02 = cos(A;x) — TT cot(7ra;) sin(A;a;), (62) 

ipi = [E - Ei{Ei)]sm{kx), (63) 

(p2 = k cos{kx — k) — IT cot(7rx) sin(A;x — k), (64) 

(^1 = [E- Ei{Ei)]sm{kx-k). (65) 

Again we have 

^(a) = iP{b) = 0, (66) 

m = m = (67) 

that makes right part of eq. ( |36l) to be equal zero that leads to evident consequences, 
similar to those, made in case (i). 

Case (iii) . The positive branch of the eigenspectrum is the following 

En = + kl kn = ^, n = 0,l,... (68) 



and negative one 



En = -y/m^ + kl kn = ^, n = 0,l,... (69) 



Eigenfunctions are 



,(n) . /, N. An) knCOS{knX) , . 

^\ = sm{knx), (t)\ = — , (70) 

En + m 

r(n) . /, N 7(n) knCOs{knX) 

= sm{knx), (t)\ = . (71) 

En + m 
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Solutions of 

hoip = Eip 

are chosen as following way 

^1 = sin(fcx), ^2 = -^!?^, = _ (72) 

+ m 

ip2[x) = sin(fcx-fc), V9i(x) = . (73) 

h — m 

And the full trace of Green function of the transformed Dirac equation is 

tr / G{x,x)dx = -— tan(fc). (74) 



Taking into account the relation (see Appendix B) 



Wfc) = g (^^,^2)%^_^, , (75) 



we get 



E ^ 2E ^/ 1 1 \ 

n=0 " n=0 ^ " 

Again the result ( 17^ is in agreement with spectral representation result. 

Constructing transformation matrix u from eigenfunctions tp^^\x) and ?/'^^-'(x) 
with the help of simple algebra one can get 

1 _ / 7rcot(7ra:/2)/2 \ 

~ 1^ -7rtan(7rx/2)/2 J' ' ^ 

~ 7r/7r\ 7r/7r\ 

ijji (x) = ijj'^ (x) — — cot ( —X ] il^i = k cos{kx) — — cot ( —x j sm{kx) , (78) 

~ TT / TT \ TT /TT \ 

ip2{x) = ^2(^) + — tan (^^^j ^^2 = kcos{kx ~ k) + — tan y^^j sin(A:x — A;). (79) 
Again we have 

MO) = Mo) = o^i>{o) = o, (80) 

<^i(l) = <^2(1) = 0^<^(1) = 0, (81) 

from which and Eq. ( l36l) it follows that transformed eigenfunctions L(f)^"-\ Lcp^^'^ 
form complete set. 

The case (iv) is similar to the case (iii). So we omit it's detail discussion. 
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5 Discussion 



We have established the relation (!35|) that connect the difference of the full traces for 
the Green functions of the initial and the Darboux transformed Dirac problems with 
energies Ei, E2 of initial states, whose wave functions are used for the construction 
of the transformation matrix and boundary values of solutions of the initial and the 
transformed Dirac equations. These relation are used to check the completeness of 
set of wave functions, obtained by the Darboux transformation of eigenf unctions of 
the initial Hamiltonian for some typical boundary problem. 
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Appendix A 

Here we will establish the relation between Wronskians of initial and Darboux trans- 
formation problem. ^ 

Using the identity 7'^7 = 1, we represent W in the form: 

W = ■?/'"^7(^ = tr[7<^'^"^] = tr[7(^'^"^7"^7] = tr[(7(^)(7'^)"^7]. (82) 

Now by usage of the definition 

ip — jp' — UxU~^il'j (83) 
(f — (f' — UxU~^(p (84) 

and the Dirac equation 

7^^' = {E-V)i;, ^<f>'^{E-V)<j>, (85) 
^yUx — uX — Vu 
we can present ■0 and <^ in the form 

(p ^ {E- fl)ip, ^={E- Q)ip, n = uXu'^. (86) 

Then we have 

W = tr[y{E -n)if^'^{E -n'^)] 

= tr[E^-f^'ilj'^ - E{'yn + n'^-f)^ilj'^ + n-fnipij'^]. (87) 

It is easy to verify that for arbitrary (2x2) matrix Q and 7 ' ^ ^ 



-1 



and 

n^7n = 7detn. (89) 

Thus, 

W = (tnip'ijj'^)(E^ - EtrQ + det Q) ^W(E^- EtrQ + det Q). (90) 

Further it is easy to check for arbitrary nonsingular 2x2 matrix u and 

Ai 
A2 



A 
that 



trfl = truAvr^ = Ai + A2, 

detO = det{uXu~^) = A1A2. 

Thus, we obtain the expression related the Wronskians of the transformed and 
the initial problems 

W = (£;- Ai)(£;- A2)W^, (91) 

where Ai = E^,, A2 = E^^- 
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Appendix B 

Here we will prove the relations (!50l) . (!75|) used in the text to check the correspon- 
dence of two results for full trace of the Green function of the Dirac problem with 
boundary conditions. 
Let us consider 



ifi(x) = r(i + x)r(i-x) 



sin(7rx) 



where T{z) is Eiler F-function. Then 
d In Hi (x) 



dx 



^p(l + x) — ^pil — x) = vr cotfvrx), 

X 



dln{r{z)) 
dz 



-etc- 't-) 

^ — ^ V'^ n — I + z J 

71=1 ^ ' 



Using the representation 

i;iz) = 

where C is the Eiler constant, it is easy to get 

1 \ ^ 2x 

77 cot nx = — h > — 
X ^-^ X — 

n=l 

By substitution x = k/n in the last relation, we have 

2k 



cot(A:) = t + J:^ 



k ' ^ k"^ — TT^n^ 

n=l 



Next we consider 



H^ix) = r(l/2 + x)r(l/2-x) 



IT 



dx 



cos(^7ra;j 

4){l/2 + x) -?/'(l/2-x) =7rtan(7rx) 



Since 



(92) 

(93) 
(94) 

(95) 



(96) 



(97) 

(98) 
(99) 



^ \n + l 

n=0 ^ 



ij{l/2 + x) - tlj{l/2 - x) 
1 



^1^ + 1 n + l/2 + xJ 

n=0 ^ ' ' 



E 



2x 



n=0 



n + 1/2)2 -a;2' 



n + 1/2- X 
by substitution x = k/n we get 

oo 



(100) 



2k 



(101) 
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